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Abstract 



In this set of lectures, we review briefly some of the recent developments 
in the study of the chaotic dynamics of nonlinear oscillators, particularly 
of damped and driven type. By taking a representative set of examples 
such as the Duffing, Bonhoeffer-van der Pol and MLC circuit oscillators, 
we briefly explain the various bifurcations and chaos phenomena associ- 
ated with these systems. We use numerical and analytical as well as ana- 
logue simulation methods to study these systems. Then we point out how 
controlling of chaotic motions can be effected by algorithmic procedures 
requiring minimal perturbations. Finally we briefly discuss how synchro- 
nization of identically evolving chaotic systems can be achieved and how 
they can be used in secure communications. 
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1 Introduction 



Dynamical systems are often modelled by differential equations when the flow is 
continuous. For an N-particle system moving under the influence of (constraint 
free) internal and external forces, Fi = 1, 2, ...N, Newton's equations of motion, 

^i^=F\r,^r,.■■■rN,^,.■.^,ty^ = l,2,...N, (1) 

where = fi(i) is the position vector of the ith particle in an inertial frame, 
subjected to 6-/V initial conditions f'i(O) and (^)(^g, constitute a deterministic 
description of the system. As long as the initial position and velocities as well 
as the nature of the forces are precisely prescribed, by solving the system of 6iV 
coupled second order ordinary differential equations (ODEs) the future evolution 
can be predicted with as much accuracy as desired. Of course the Laplace dictum 
of complete determinism fails when statistical fluctuations and quantum effects 
are present. Barring the presence of these forces, one might expect that the 
solution of the initial value problem of (1) should lead to complete predictability. 
However, in recent times it has been clearly demonstrated that a new kind of 
indeterminism or unpredictability can arise when the forces acting on the system 
are nonlinear (even when statistical and quantum forces are absent) in suitable 
form. 

When the forces Fj in (1) are linear, we only have to solve a system of linear 
differential equations wherein the superposition principle holds good. Any small 
deviation in the initial conditions or numerical errors will then grow at the 
worst only linearly, so the dynamics of the system can be exactly predicted. 
However when the forces Fi are nonlinear, one has to solve nonlinear differential 
equations. Then for suitable nonlinearities, one typically finds the possibility 
that any small deviation in the prescription of initial conditions or round-off error 
at any stage of the calculation can multiply quickly, leading to an exponential 
divergence of trajectories corresponding to nearby initial conditions. Or, in other 
words, the system can be extremely sensitively dependent on initial conditions 
leading to complex or chaotic motions. Physically this means that any small 
fluctuation during the evolution can lead to realizable macroscopic effects, the 
so-called 'butterfly effect' [1,2]. Mathematically, the phase trajectory might 
end up asymptotically in a strange attractor of zero volume but of fractional 
dimensionality in the case of dissipative systems, or it can fill the whole phase 
space densely in the case of conservative systems. 

The notion of complexity and chaos can be realized even in low-dimensional 

nonlinear systems. Particularly, damped and driven oscillators are physically 
realizable examples exhibiting the immense varieties of bifurcations and chaos. 
They are often modelled by a single, second-order, nonlinear differential equation 
(or equivalent first order system) with periodic inhomogeneities. These include 
the Duffing oscillator, the Bonhoeffer-van der Pol oscillator, the Duffing-van 
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der Pol oscillator, the Murali-Lakshmanan-Chua electronic circuit, the damped 
driven pendulum, the driven Morse oscillator and so on [3,4]. In this set of 
lectures, wc will consider the extremely varied bifurcations exhibited by these 
systems and the various types of regular and chaotic motions admitted by them 
and some of their applications. 



2 Linear and nonlinear oscillators 



There are obvious characteristic differences between the dynamics of linear and 
nonlinear oscillators. We will briefly consider these aspects with simple exam- 
ples. 



2.1 Linectr oscillators 

We will point out how predictability becomes an inherent property of linear 
oscillator systems. We will also consider the type of attr actors exhibited by 
these systems in this section. 

2.1.1 Free linecir harmonic oscillator 

The ubiquitous linear harmonic oscillator of unit mass satisfying the equation 
of motion ^ 

^+u;lx = 0, u;l>0 (2) 
and the initial conditions a;(0) = A, {dx/dt)t=Q = admits the solution 

x{t) = Acosuot (3) 

such that 

E=lx^ + 1„„V = 1.SA' (. = i) (4) 

is a constant. The corresponding phase trajectories in the phase space {{x-x)- 
space) are concentric ellipses with the origin (0, 0) being a centre type (elliptic) 
fixed point (Fig. la). 



2.1.2 Damped linear harmonic oscillator 

When viscous damping force is present, the harmonic oscillator equation becomes 

— +a-+.o^. = 0, (5) 
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where a (> 0) is the damping parameter. For the same initial conditions, a;(0) = 
a;(0) = 0, the solution for the physically interesting underdamped case (a < 
2a;o) becomes 

where 

Solution (5) represents a damped oscillatory motion and energy is dissipated. 
The corresponding phase trajectory is a spiral asymptotically approaching the 
fixed point at the origin which is now a stable focus (Fig. lb). 



2.1.3 Damped and driven linear oscillator 

When an external sinusoidal force is added to the system, we have 

(Px dx o , . 

Again for the same set of initial conditions as above, we have the solution 



VaV4 + C2 f-at 

where 



x{t) = At— — exp [ j cos{Ct — 5) + Ap cos{(jt — 7), (7a) 



f -1 / ttw \ 

and At is chosen so as to satisfy the initial condition x{0) = A, x(0) = 0. The first 
term in (7) is obviously a transient which decays for large t. Thus asymptotically 
the system oscillates with a frequency equal to that of the external periodic 
force. At the resonance value co = ivo, the amplitude of oscillation Ap takes the 
maximum value. The associated phase trajectory corresponds asymptotically 
to an attractor which is now a limit cycle as depicted in Fig. 1(c). In all 
the above cases, the initial value problem is exactly solvable, the motion is 
completely regular and the future is fully predictable, given the initial state. 
Trajectories corresponding to nearby initial conditions continue to be nearby, 
leading to complete predictability. The attractors are only fixed points or closed 
curves corresponding to regular motions. 



2.2 Nonlinecir oscillators 



As a prelude to understanding the typical behaviour of damped and driven 
nonlinear oscillators, we ask how the dynamics of the previously studied linear 
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oscillators are modified when a nonlinear spring force of cubic form is added so 
that the equation of motion becomes 



^ +a'^ +colx + Px^ = fsmojt. (8) 



2.2.1 Forced oscillations - Resonances (/ ^ 0) 

One immediately observes that the behaviour of nonlinear oscillators of the form 
(8) can be qualitatively quite different from that of the linear case, even for very 
small (3. Near the primary resonance oj Ri uiq, /3 <C 1, one can write the lowest 
order approximate solution in the form 

x{t) = Acos{ut + 5), (9) 

where ^ is a constant to be fixed. Using (9) in (8) one obtains the resonance 
curve 

+ {aAcof = f. (10) 

A typical form of this curve is shown in Fig. 2 for f3 > 0, which clearly shows 
the associated hysteresis phenomenon. In addition one also observes secondary 
resonances (subharmonic and superharmonic) at frequencies w w nivo, and u) = 
^,n = l,2,3,... 



2.2.2 Nonlinear oscillations and bifurcations [1-7] 

However, to obtain a complete picture of the dynamics, one has to have recourse 
to detailed numerical analysis, using for example the fourth order Runge-Kutta 
integration method. As a result one obtains a rich variety of bifurcation phe- 
nomena, namely successive qualitative (and quantitative) changes in the nature 
of oscillations as the nonlinearity parameter (or another control parameter say / 
or Cij or a for fixed nonlinearity parameter) is varied. A mind-boggling multitude 
of bifurcations leading to the onset of chaos and further bifurcations occurs in 
the parameter space. The details are further described in Sec. 4. The above 
type of numerical analysis can be used in different ways to extract information 
about the dynamics. Some of them are the following. 

1. Trajectory plot ((.T-t)-plot) 

2. Phase portrait: (a;-i;)-projection of the trajectories from three-dimen- 
sional {x,x,t)-space. 

3. Poincare map: Stroboscopic (snapshot) portrait of the trajectory at 
every period T(= 
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4. Power spectrum: P{u}) 



^ I x{t)exp{-ii2)t) 



5. Bifurcation diagram: Plot of Poincare points of the solution x{t) vs 
control parameter. 

6. Lyapunov spectrum: Plot of maximal Lyapunov exponent vs control 
parameter. 



Besides these, one often uses correlation functions, basins of attractors and so 
on. 

What are the kinds of behaviours a dissipative nonlinear system of driven 
type typically exhibits and how do they occur? The common types of behaviours 
are the following: 

1. Fixed points 

2. Periodic response of period w of the input 

3. Subharmonic response (fi = ^, n integer) 

4. Superharmonic response {Q — nw) 

5. Almost periodic (quasiperiodic) response (consisting of periodic com- 
ponents whose frequencies are incommensurable) 

6. Chaotic response (nonperiodic and sensitively dependent on initial 
conditions) 



Some typical attractors are illustrated in Figs. 3. There are typical routes 

through which periodic and chaotic responses can arise as the control parameter 
varies. Three most predominant routes arc the following: 

1. Fcigenbaum's period doubling scenario. 

2. RucUe-Takcns-Nowhousc scenario of a quasiperiodic route. 

3. Pomeau-Manneville scenario of an intermittency route. 



Besides these three prominent routes, there are many other less common routes, 
such as period adding sequence, equal-periodic bifurcations and so on, which 
occur occasionally in dissipative systems. 



2.3 Linear and nonlinear electronic circuits [6,8] 

Besides numerical and analytic studies of nonlinear oscillators there is another 

useful way of investigating them, namely by analog simulation through suitable 
electronic circuits. Such studies help in a dramatic way to scan the parameter 
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space quickly and to avoid long transients as in the case of numerical studies. 
The construction of such circuits is fairly easy in typical cases, giving rise to 
quick responses as the control parameter is varied. 

Prom another point of view, a variety of nonlinear electronic circuits consist- 
ing of nonlinear physical devices such as nonlinear diodes, capacitors, inductors 
and so on or devices constructed with piecewise- linear circuit elements have been 
utilized in recent times as model nonlinear systems in their own right to explore 
different aspects of chaotic dynamics. They further give rise to interesting tech- 
nological applications. 



2.3.1 Linear circuits 



Linear circuits consist of only linear elements such as linear resistors, capacitors 
or inductors besides the source. In fact, the ubiquitous LCR circuit (Fig. 4) 
can be considered as a prototype analog simulation of the damped and driven 
linear oscillator (8). The LCR circuit of Fig. 4 consists of a linear resistor, a 
linear inductor, a linear capacitor and a time-dependent voltage source, f{t) = 
Fg sinwgt. Applying Kirchhoff's voltage law to the circuit of Fig. 4, we obtain 

L— + RiL + v = Fssinwst, (11) 

with ii,(0) = io, v{0) = vq, where Zl is the current along the inductor, v is 
the voltage across the capacitor, L is the linear inductance and R is the linear 
resistor. 

Now substituting = into Eq. (11) and rearranging one obtains the 
linear second-order inhomogeneous ODE with constant coefficients, 

(Pv Rdv 1 Fs . 

Obviously Eq. (12) is analogous to (6) with appropriate identifications, and its 
solution can be immediately written down in the form (7) corresponding to limit 
cycle oscillation with frequency co = LOg- The actual experimental result for the 
wave form v{t) and the phase portrait in the (t;-ii)-plane is shown in Fig. 5, 
which is seen to mimic the solution of (6) given in Fig. 1(c). 



2.3.2 Nonlinear circuits [9] 
a) analog simulation 

As mentioned above, nonlinear oscillator equations can be analog simulated by 

the use of suitable operational amplifier modules and multipliers. Normally 
op-amps are used as integrators, differentiators, sign-changers, adders and so 
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on. The mathematical operations are generally obtained cither by using an op- 
amp individually with resistor feedback, or an op-amp with capacitor feedback. 
Similarly, product and division operations can be performed by multiplier and 
divider chips such as AD532 or AD534. One such analog simulation circuit 
constructed with adding integrators, scale changers and multipliers allows us 
to study the chaotic dynamics of the cubic Duffing oscillator equation (8). Its 
schematic diagram is given in Fig. 6. 

b) Piecewise-linear circuits: MLC circuit [10] 

Considering the familiar LCR circuit given in Fig. 4, let iis add a nonlinear 
resistor, the so called Chua's diode characterized by a three segment negative 
resistance characteristic curve as given in Fig. 7. The modified circuit is given 
in Fig. 8. This simple non- autonomous circuit was first reported by Murali, 
Lakshmanan and Chua [10]. The equations of motion governing this circuit are 

C-=iL-g{v) (13a) 

L-^ = —V — RiL + FgSmuJst, (136) 

where g{v) = GbV + \{Ga-Gb)[\v + Bp\ — \v — Bp\], and Ga, Gb and Sp are circuit 
parameters. It can also be considered as a prototype of the chaotic nonlinear 
oscillators. 



3 Some ubiquitous damped and driven nonlin- 
ear oscillators [6] 

3.1 Duffing oscillator 

One of the widely studied damped and driven nonlinear oscillators is the Dufiing 
oscillator given by Eq. (8) . It can be considered to represent a particle (of unit 
mass) moving in a potential well 

V = lu;lx' + (15) 

subject to viscous damping and external periodic force. Here one may distinguish 
three cases. 

1. Double- well (twin- well) potential loq < 0,/3 > 

2. Single-well potential ljq > 0,l3 > 

3. Double-hump potential: > 0, /3 < 0. 
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Each one of the above three cases is an important model for many physical 
problems. The double-well potential represents the oscillations of a buckled 
beam under a periodic force, plasma oscillations and so on. Similarly, the single- 
well oscillator is a standard model of a particle oscillating under viscous damping 
and a periodic force, while a double-hump potential well is a model for weakly 
pinned charge density waves in anisotropic solids and superionic conductors. 

Now Eq. (8) has an obvious sealing property which allows one to fix two of 

the five parameters a, ciJ§, /?, / and oj. From an experimental point of view it is 
more convenient to fix the spring forces Wq and (3 and vary the parameters in the 
(a-a;-/)-space. Prom a practical point of view one normally studies behaviour 
in the /-space or (/-a')-space by fixing the remaining parameters as well. In the 
next section (Sec. 4), we will briefly discuss the chaotic dynamics of all the three 
potential wells of Eq. (8). 



3.2 BVP and DVP oscillators 



The nerve impulse propagation along the axon of a neuron can be studied by 
combining the equations for an excitable membrane with the differential equa- 
tions for an electrical core conductor cable, assuming the axon to be an infinitely 
long cylinder. A well known approximation of FitzHugh and Nagumo to describe 
the propagation of voltage pulses V{x, t) along the membranes of nerve cells is 
the set of coupled PDEs 

-Vt = F{V) + R-I, (16a) 

Rt ^ c{V + a~bR), (16&) 

where R{x, t) is the recovery variable, / the external stimulus and a, b, c are re- 
lated to the membrane radius, specific resistivity of the fluid inside the membrane 
and temperature factor respectively. Here suffixes stand for partial differentia- 
tion. 



When the spatial variation of V, namely Vxx, is negligible, Eq. (16) reduces 
to a set of two coupled first-order ODEs known as the Bonhoeffer-van der Pol 
oscillator system, 

V = V- — -R + I (17a) 

o 

R=c{V + a-bR), (—^^ (17&) 

with F{V) = — V + Normally the constants in Eq. (16) satisfy the inequal- 
ities 6 < 1 and 3a + 2b > 3, though from a purely mathematical point of view 
this need not be insisted upon. Then with a periodic (ac) applied membrane 
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current Ai coswt and a (dc) bias Aq, the BVP oscillator equation becomes 

V = V - — - R + Ao + AicosLut, (18a) 
o 

R = c{V + a^bR). (186) 

Again system (18) shows a very rich variey of bifurcations and chaos phenomena 
as discussed in Sec. 5 below. 

Further, Eqs. (18) can be rewritten as a single second-order differential 
equation by differentiating (18a) with respect to time and using (18b) for R, 

( ^ ■ be 

V-(l-bc) \ 1 } V-c(b-l)V+ — V^ = c(Aob-a)+Ai cos(u)t+(p), (19) 

[ 1 — oc J 3 

where (f) = tan"^ ^. Using the transformation x = {1 — bc)~^^^'^^V, t — > t' = 
t+ ^, Eq. (19) can be rewritten as 

X + p{x^ - l)x + ujIx + I3x^ = fo + ficos ut, (20) 



where 



p=(l-6c), w2 = c(l-6), 
P = bc , /o = c- 



3 ' ' 

Eq. (20) or its rescaled form 

x + p{kx^ +g)x + ujIx + I3x^ = fQ + fi cos uit (22) 

is the Duffing-van der Pol equation. In the limit k=0, we have the Duffing 
equation discussed above (with /o = 0), and in the case (3 — (g — ^1, k = 1) 
we have the driven van der Pol equation. Eq. (22) again exhibits a very rich 
variety of bifurcations and chaos phenomena, including quasiperiodicity, phase 
lockings and so on, depending on whether the potential V = \u>qX^ + is i) a 
double well, ii) a single well or iii) a double hump. We will discuss some of the 
details in Sec. 5. 



4 Duffing oscillator: Bifurcations and chaos 

As noted earlier, the Duffing oscillator (Eq. (8)) dynamics can be described 
by a phase diagram in a three-parameter (a-/-a;)-space. In the following, we 
will concentrate on the /- or w- or (/-a;)-parameter space because of practical 
considerations. 
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4.1 The double-well oscillator {ujI < 0,/? > 0) 



Considering a particle in a double-well potential V = —1/2\uj'q\x'^ + with a 

base vibrating periodically, wc notice that there are three equlibrinm points in 
the force-free case (/ = 0), corresponding to Lo^x + fix^ = 0. They are the stable 
fixed points (centres) 

= ±yS (23) 

and an unstable fixed point which is a saddle, 

a^o = 0. (24) 

Then the oscillations about the two stable equilibrium points are given (after a 
redefinition of x) by 



x + ax + 2\ujl\xT^\l\^l\li x^ + I3x^ = f sin Lot. (25) 
4.1.1 Period doubling scenario in /-space 

Choosing the initial condition such that the system approaches the left equilib- 
rium point {xx = —1.0) asymptotically as example, and choosing the parameters 
as (jJq = —1.0, (} = 1.0, a = 0.5, to = 1.0, we can solve Eq. (25) numerically 
and obtain the following picture. To start with we realize from the theory of 
nonlinear oscillations that the system can exhibit two types of periodic motions: 

1. Small orbits {SOs): Oscillations around the stable fixed point xi. 

2. Large orbits (LOs): Large amplitude oscillations that encircle all the 
three equlibrium points xi, X2 and xq. 

Then the following picture arises. 
A. Small orbit oscillations 

As / is increased slightly from zero, a stable period T limit cycle about the 
left equilibrium point x\ = — 1.0 occurs, which persists up to / = 0.34. Then a 
period doubling scenerio starts at / = 0.35 with a 2T periodic cycle, followed by 
a AT periodic cycle at / = 0.358, etc. This cascade of bifurcations accumulates 
at / = /c = 0.361, the sequence converging geometrically at the rate known as 
the Feigenbaum ratio. Beyond / = /c, one finds chaotic orbits occurring where 
the motion is highly dependent on the initial conditions. This can be checked by 
the various characterizations discussed previously. These motions are illustrated 
in Figs. (9a-c). 
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B. Hopping Oscillations 

When / = 0.37, the onset of a crisis occurs, in which the chaotic motion 
that had been confined to the left weU suddenly expands to cover the right well 
also, in a hopping cross-well chaotic motion. 

C. Large orbit oscillations 

For still larger values of /, the chaotic state suddenly disappears and is 
replaced by a stable period-mT oscillation, overcoming the central barrier, ex- 
tending over both the valleys and encompassing all the equilibrium points. (Fig. 
9f, g). 

Fuller details are given in Table 1 and the bifurcation diagram in Fig. 10. 

D. Antimonotonicity 

In addition to the periodic window resulting from the crisis, there is also 
another possibility due to antimonotonicity. Here one finds that periodic orbits 
are not only created through period doubling cascades but are also destroyed 
through reverse period doubling when a control parameter is monotonically in- 
creased in any neighborhood of a homoclinic tangcncy value. Two types of such 
antimonotonicity are shown in Fig. 11. Fig. 11(a) represents a period-3 bubble 
in the range 0.53 < / < 0.62, while Fig. 11(b) shows a reverse-period bubble in 
the region 0.43 < / < 0.48. 

4.1.2 Phase-diagram in (/-a;)-space 

To appreciate the rich variety of bifurcations, let us look at the (/-w)-phase 
diagram (Fig. 12) given by Szemplinska-Stupnicka and Rudowski [11] in the 
frequency range 0.25 < uj < 1.1, fixing a = 0.1, |t^oP = 0.5 and (3 = 0.5. 
One observes that SO, symmetric LO, unsymmetric LO and cross- well chaotic 
or stable attractors coexist. The SO motions occur within the whole [f-uj)- 
planc except for two T/-shapcd regions, one with a cusp at w = 0.8 and / = F2 
at the principal resonance region, the other cusp at w = 0.4 , and / = 0.14, 
that is the superharmonic zone. Inside the two F-shaped regions the system 
can exhibit cross-well chaotic (or regular) motion. Also, in some regions, two 
different attractors coexist while, in the others, a single steady state (globally 
stable) motion can be observed. Typical orbits: the form of the various single 
and coexisting states denoted in Fig. 12 by numbers 1, 2, .... 20 are shown in 
Fig. 13. Regular attractors are illustrated by their phase portraits and chaotic 
attractors by Poincare maps. They include all the types of motions discussed 
above. 
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4.2 Single well and double hump potentials 



Considering the single-well potential, one observes hysteresis in the primary 
resonance region, harmonic resonances and symmetry breaking orbits for a range 
of low values of /. For sufficiently high values of /, chaos occurs through a 
period-doubling route. Similarly, in the double-hump case, one finds that one 
symmetric attractor bifurcates into two mutually symmetric attractors. Both 
attractors then undergo period-doubling cascades to chaos, where windows of 
both even and odd periods occur. Then both chaotic attractors broaden and 
merge to form a single attractor, ultimately leading to a boundary crisis. For 
further details see Ref. [6]. 

4.3 analog simulation and experimental verification 

The Duffing oscillator (8) can be easily analog-simulated using op-amps and 

four-quadrant multipliers. Replacing x(t) = J xdt, x{t) = J xdt, the equivalent 
schematic simulation circuit is given in Fig. 6. Here blocks ICi and IC2 rep- 
resent two integrators, ICs is an inverter (sign-changer), and Mi and M2 are 
four-quadrant multipliers. The experimental results for the double- well case are 
given in Figs. 14, and they agree with the theoretical numerical analysis. For 
further details see [6] again. 



5 Dynamics of the BVP and DVP oscillators 
[6] and MLC Circuit 



Let us now quickly review the chaotic dynamics and bifurcations of the BVP 
and DVP oscillators discussed in Sec. 3.2. Apart from the period-doubling route 
to chaos, one also observes the presence of phase-lockings and devil's staircase 
structures, quasiperiodic and intermittency routes to chaos predominantly. The 
nonlinear dynamics of the MLC circuit are also discussed. We shall therefore 
give only a very brief account of these structures. 

5.1 BVP oscillator 

Considering Eq. (18), the equilibrium points in the absence of external stimuli 
{Ao = Ai = 0) are the roots of the equations 



V^-qV-p = 0, 



(26a) 



V + a-bR = 0, 



(266) 
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where p = and q = 3-^^^^^. Eq. (26a) has three real roots for 27p^ > and 
one real root for 27p'^ < q^ . The linear stability analysis of Eqs. (26) shows that 
for a fixed point {Vq, Rq) to be stable, the necessary and sufficient condition is 

{Vi + 6c - 1) > 

and 

c(l + hV^ - 6) > 0. (27a) 

Thus for suitable choices of parameters one can have one stable focus or one 
saddle and two stable foci. 

Proceeding along the same lines, with the inclusion of constant bias {Ai = 0, 
Aq > 0), one can find fixed points by simply modifying p to p = 2>{Aq — a/h) in 
Eqs. (26). It is of particular interest to check whether the fixed points undergo 
a Hopf bifurcation. This can be easily verified by checking the nature of changes 
in the stability property of the fixed points and the parametric choices for which 
the eigenvalues cross the imaginary axis. For the present problem, one can easily 
check that a Hopf bifurcation occurs at 

= % - + I {27b) 

One checks that for < ^o- and Aq > Aq+ the fixed points are asymptot- 
ically stable, being either a node or a spiral point. For Aq G (Ao_,Ao+), the 
eigenvalues have positive real parts so that the fixed point is unstable and one 
may expect a stable limit cycle in this range. 

Similarly when only a periodic input current {Ai 0, Aq = 0) is present, 
the system, as in the case of the Duffing oscillator, undergoes period-doubling 
bifurcations to chaos, followed by windows, intcrmittcncy and so on. The bi- 
furcation diagram given in Fig. 15 describes briefiy the various bifurcations 
admitted by this system for the parametric choice a = 0.7, b = 0.8, and c = 0.1. 

Now if we consider the situation in which both external d.c. and periodic 
stimuli are present, the system enters into an interesting phenomenon called 
phase-locking or mode- locking. Here the natural frequency of the oscillator (un) 
and the frequency of the driving force get locked. When one calculates the 
winding number W = ^ &s a, function of the d.c. current Aq (which will give 
rise to a different f„), the corresponding schematic diagram has a devil's staircase 
structure. A typical example for the choice Ai = 0.74 is shown in Fig. 16. 

5.2 DVP oscillator 

One can consider three different cases for the DVP oscillator (22), as in the case 
of the Duffing oscillator, namely i) double well, (ii) single well and (iii) double 
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hump. They exhibit many comphcated bifurcation phenomena including period- 
doubhng, hysteresis, quasiperiodicity, phase-locking and so on. Some interesting 
features of the double-hump DVO oscillator include 

i) the existence of a Farey sequence, 

ii) double hysteresis loops, 

iii) a swallow-tailed bifurcation structure of subharmonic locking, 
and so on. For more details see Ref. [6]. 

5.3 Chaotic dynamics of the MLC circuit [10] 

The simple second-order nonautonomous dissipative nonlinear circuit consisting 
of Chua's diode as its only nonlinear element, suggested by Murali, Lakshmanan 
and Chua (MLC) , can exhibit a rich variety of bifurcations and chaos phenom- 
ena. The circuit diagram and its state equations have already been introduced 
in Sec. 2.3.2. Rescaling Eq. (13) with v = xBp, il = G.y.Bp, G = 1/R, cj = ^ 
and t = TCG, and then redefining T as t, we have the following set of normalized 
equations, 

x = y — h{x), y = —l3y — ul3y — l3x + f sinwt, (28a) 
where /? = u = GR, and / = Here h{x) stands for 

{bx + a — b, X > 1 
ax, \x\ < 1 (286) 
bx — a + b, X < —1. 

Here a = 6 = GbG. Now the dynamics of Eq. (28) depend on the parameters 
p, (3, a, 6, to and /. The experimental circuit parameters used in the actual 
circuit correspond to P = 1, v = 0.015, a = —1.02, 6 = —0.55 and to = 0.75. It 
is straightforward to establish that there exists a unique equilibrium {xq, yo) for 
Eq. (28) in each of the following three regions, 

D, = {{x,y)\x>l}: P+ = (-fci,-fc2) 
Do = {{x,y)\\x\<l}: O = (0,0) 
D_i = {{x, y)\x<-l}: P" = (fci, ^2) 

where h = [((T(a - b))/{p + ab)], fcs = [/3(6 - a)/(/3 + ab)] and ct = /3(1 -h u). 
In each of the regions Di, Dq, D_i, Eq. (28) is linear. One can easily check 
by a stability analysis that O G Dq is a. hyperbolic (saddle) fixed point, while 
and P~ are stable spiral fixed points. Depending on the initial conditions, 
any of them can be observed. As the forcing signal / is increased (/ > 0), the 
stable fixed points undergo a Hopf bifurcation. A further increase in / leads to 
period doubling bifurcations and finally to chaos. The corresponding numerical 
and experimental results are given in Fig. 17. 
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6 Dictionary of Bifurcation 



The various oscillator systems discussed in the previous sections exhibited a vari- 
ety of oscillatory behaviours including fixed points, resonant periodic oscillations, 
period-doubled oscillations, quasiperiodic and chaotic motions. In particular, the 
various bifurcation diagrams clearly show that all these systems exhibit quali- 
tative (and quantitative) changes at critical values of the control parameters as 
they are varied smoothly. These correspond to the different bifurcations exhib- 
ited by the systems. The question is how to understand the existence of such 
bifurcations and the mechanism by which they arise. One possible approach is 
to look for the local stability properties of solutions in the neighbourhood of the 
critical parameter values at which bifurcations occur. However, since no exact 
solutions are available for the above systems, one can prepare a dictionary of 
bifurcations [12] that occur in simple and interesting low dimensional nonlinear 
systems, and then use them as reference systems for the identification of the 
bifurcations in the various nonlinear oscillators. In Tables 2 and 3 we give their 
salient features which are self-explanatory. 



7 Analytic approaches 



The various types of bifurcations exhibited by a typical nonlinear oscillator can 
be classified by using the dictionary of bifurcation phenomena like the ones dis- 
cussed in the earlier section. However, some approximate analysis can be per- 
formed analytically by constructing pcrturbativc solutions and studying their 
stability properties. Moreover the analytic structure of solutions in the com- 
plex f-plane can throw much light on the integrability and nonintegrability and 
chaotic behaviour of dynamical systems. We will discuss some features of these 
aspects in this section. 



7.1 Theoretical analysis of large orbit T-periodic solution 
in the double-well Duffing oscillator 

The numerical analysis presented in Sec. 4 for the double-well Duffing oscillator 
conclusively showed that it admits a large orbit (LO) T-periodic solution for a 
wide range of drive frequencies. One might then expect a first-order approximate 
solution whose stability analysis might give a good estimation of the paranic;ter 
domain in which it occurs. So we shall look for a T-periodic solution of the 
oscillator equation [11], 

—^+a—-\ujl\x + (ix^=fcoswt, (29) 
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which is close to a harmonic function of time as 

a;(0) = A{0) cos(wt + (30) 

where ^(0) and </> are to be determined. For this purpose we transform (29) into 
the form 

x + oj'^x + ef{x,x,t) =0 (31) 

as 

X + uP'x + \x{ax + xa — fl^x — \u)o\'^x + j3x^ — f cos tot) = 0, (32) 

where 

fia = a, jjLu = jjtVt^ - , l^f = f, mI'^oI^ = |wo|^,M/^ = /3- (33) 

Here the quantity fi has been introduced from the form of the harmonic solution 
X = acosV,t, = —\u)q \ + 3/4/30^ > 0. Now one can look for a periodic solution 
of (32) in the form 

x{t,fi) = xo(ro,ri, . . .) + nxi{To,Ti, . . .), (34) 

oo 

where the multiple scales To = t, Ti = fit, . . . and ^ = jJ'^Dn, Dn = so 

n=o " 

that we have the system of equations 

DqXq + lo^xq = 0, (35a) 

DqX\+uj'^X\ = —2DoDxXQ — aDoXo — {a — ^^—\uiQ)\^)xo — (}xQ + fcosujt, (356) 

etc. The general solution of (35a) is 

Xa = Ai(Ti) cxpio^To + C.C., (36) 

where Ax is, in general, complex. Using (36) in (35b), in order to eliminate the 
secular terms in the solution of X\ , one requires that 

-[2iLo{DiAi + a/2Ai) + {a-if - |wo|^)^i + 3l3AlAl] + f/2 = 0. (37) 

Making the substitution 

Ai(Ti) = (l/2)a(Ti) exp#(Ti) (38) 

and equating real and imaginary parts for steady state conditions with w 

+ 3/4/3a2, we obtain 

Finally, the first-order correction becomes 

txxi{t)=A3Cos3{iJt + <j>), ^3 = 1^ = ^, (40) 
so that the LO solution becomes 



x{t) = acos{ujt + </>) + -— -o cos3(a;t + (j>) = x{t + T). (41) 
32cj 
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7.1.1 Stability of LO solution, 



A. Soft mode instability 

Considering small disturbances to the amplitude and phase of the solution 
(41) leading to 

x{t) = (a + 5a) cos(wt + 4) + 5(1))+ ^"g^^"^ cos 3(wt + 4) + 5(1)), (42) 
one obtains the linear eigenvalue problem 

-f )W)- 

Then the linear stability of the solution (41) depends upon the (2 x 2) matrix 
in the right hand side of (43) from which a necessary criterion for the stability 
of the solution can be obtained (see Fig. 20). 

B. Symmetry-breaking instability 

Considering more general instabilities, x{t) = x{t) + 5x{t), one can deduce 
a Hill type eigenvalue problem, 

5x + aSx + {Xo+ ^ Xn cos n(jt)5x = 0, (44) 

n=2,4,6 

where Aq, A2, A4, Ae are appropriate coefficients. One can check that while there 
is no period-doubling instability, symmetry-breaking instability can occur here. 
This may be verified by assuming a solution of the form 

5x{t) = e^i*(6o + 621 cos2a;f + 622 sin2a;t), (45) 

and finding the values of ei. The corresponding curve cossia) is given in Fig. 20. 
One finds that the curves t<;s(a) and wssio) determine the region of existence 
of LO to a reasonable approximation. 



7.2 Stability analysis of SO motion 

The numerical analysis of the double-well DufSng oscillator discussed earlier 
clearly shows the presence of SO periodic motion confined to one of the wells. 
One can therefore again look for T- periodic solution in the right well for example 
and apply multiple scale perturbation theory. Comparing with Eq. (29), we 
consider the equation 

x + ijp'x + \mjJ\x^ + [^(ax + 0)30;^ + ax — f coscjt) = 0, 
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where 



IMVi = 3/2, ii^u)3 = 1/2, ii^a = a, a\j^ = 1 — ijp' , 



(46) 



Proceeding as in the previous case one obtains the second-order approximate 
solution 

x{t) = acos{uit + (/))- 3/4a2 + l/4a^ cos2(a;t + (j)) + a^/16cos3(a;t + 0). (47) 

Considering, as in the case of LO periodic solutions, the stability of the SO 
periodic solution (47) , one can analyse the first-order instability as well as period- 
doubling instability. As a result one can conclude that the approximate T- 
periodic solution is unstable within a region of the driving frequency uia < < 
u)PD and that no other SO stable solution is possible. The above result is found 
to give rise to a surprisingly good estimation of the system's critical parameter 
values for which cross- well chaos really occurs. For more details, see Ref. [11]. 

7.3 Analytic structure of the Duffing oscillator 

An altogether different approach to investigating the nonintcgrability aspects of 
the DufHng oscillator and other dynamical systems is the study of the analytic 
structure of singularities of the solutions in the complex time plane, known 
in the literature as the Painleve singularity structure analysis [see for example 
Refs. 13-15]. Originally this method was applied successfully to isolate the 
integrable cases of the rigid body motion in a gravitational field. Its recent 
renaissance started with the discovery of soliton systems (see the lectures of Prof. 
Ablowitz in this School). One essentially tries to isolate the parameters of the 
nonlinear differential equations underlying the dynamical system for which the 
solution is free from movable critical points (movable branch points and essential 
singularities) so that the solution is meromorphic in the complex time plane. 
For nonintegrable systems, the singularities are often located in a complicated 
pattern, giving an indication that a sufficient number of analytic integrals is 
lacking. Consider for example the unforced, undamped (free) Duffing oscillator, 



The solution to (48) can be given in terms of Jacobian elliptic functions. 



where the modulus k = (iA^ /2{uj'^ + (iA^). One knows that the elliptic func- 
tions are meromorphic and doubly periodic. The only singularities they possess 
in the complex f-plane are poles of order 1 at the lattice points. 



X + ujIx + px^ = 0. 



(48) 



X 



(t) = Acn{nt + 6),n= + f3A^ 



(49) 



(Oi + 5) = 2nK + (2m -|- l)iK' , n,meZ. 



(50) 
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We also note that despite the existence of these poles, the system possesses the 
analytic integral, E = ^uJ^A^ + jf3A^, in the complex i-plane and the system is 
intcgrablc. A typical pattern of the lattice of singularities of the solution of (48) 
is given in Fig. 21a. The situation becomes more complicated when damping as 
well as forcing are present. Typical singularity structures are shown for the case 
when damping alone is present in Fig. 21b, and when forcing is also present in 
Fig. 22. One observes a progressive clustering of singularities off the real axis 
as the forcing strength increases. Actually, one can identify much more intricate 
structures through local analysis or numerical analysis. Around each singular 
point there is an infinite-sheeted, four-armed structure of singularities. For more 
details, see Refs. [13-15]. Similarly, infinite-sheeted structures are present in the 
case of the DVP oscillator [16]. 



8 Controlling Chaos 



Eventhough the presence of chaotic behaviour is generic and robust for suit- 
able nonlinearities, ranges of parameters and external forces, there are practical 
situations where one wishes to avoid or control chaos so as to improve the perfor- 
mance of the dynamical system. Also, eventhough chaos is sometimes useful as 
in a mixing process or in heat transfer, it is often unwanted or undesirable. For 
example, increased drag in flow systems, erratic fibrillations of heart beats, ex- 
treme weather patterns and complicated circuit oscillations are situations where 
chaos is harmful. Clearly, the ability to control chaos, that is to convert chaotic 
oscillations into desired regular ones with a periodic time dependence would be 
beneficial in working with a particular system. The possibility of purposeful 
selection and stabilization of particular orbits in a normally chaotic system, us- 
ing minimal, predetermined efforts, provides a unique opportunity to maximize 
the output of a dynamical system. It is thus of great practical importance to 
develop suitable control methods and to analyze their efficacy. Very recently 
much interest has been focussed on this type of problems [6,17-19]. 



8.1 Algorithms for controlling chaos: Feedback and non- 
feedback methods 

Let us consider a general n-dimensional nonlinear dynamical system, 

X = ^=F{X ■ p; t), (51) 

where X = {xi,X2,X3, ...,a;„) represents the n state variables and p is a control 
or external parameter. Let X{t) be a chaotic solution of Eq. (51). Different 

control algorithms are essentially based on the fact that one would like to effect 
the most minimal changes to the original system so that it will not be grossly 
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deformed. From this point of view, controlling methods or algorithms can be 
broadly classified into two categories: 

(i) feedback methods, and 

(ii) non- feedback algorithms. 

Feedback methods essentially make use of the; intrinsic properties of chaotic 
systems, including their sensitivity to initial conditions, to stabilize orbits al- 
ready existing in the systems. Some of the prominent methods are the following. 

1. Adaptive control algorithm 

2. Ott-Grebogi-Yorke (OGY) method of stabilizing unstable periodic 
orbits 

3. Singer et al.'s method of stabilizing unstable periodic orbits 

4. Various control engineering approaches. 

In contrast to feedback control techniques, non-feedback methods make use 
of a small perturbing external force such as a small driving force, a small noise 
term, a small constant bias or a weak modulation to some system parameter. 
These methods modify the underlying chaotic dynamical system weakly so that 
stable solutions appear. Some of the important controlling methods of this type 
are the following. 

1. Parametric perturbation method 

2. Addition of a weak periodic signal, constant bias or noise 

3. Entrainment-open loop control 

4. Oscillator absorber method. 



For more complete details of the various controlling methods, see for exam- 
ple Refs. [6,17,18]. 

Here is a typical example of adaptive control algorithm. We can control the 
chaotic orbit Xs = {xsjUs) of the EVP oscillator by introducing the following 
dynamics on the parameter Ai: 

X = X — - — y + Ao + Ai cos uit, (52a) 

y = c{x + a - by), (526) 

ii = -e[{x - Xs) - (y - ys)], e « 1 (52c) 

The actual conversion of the chaotic motion into period-T and 2T limit cycles 
in a typical case is shown in Fig. 23. All the other methods can be discussed in 
similar fashion. 
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9 Synchronized Chaotic Systems and Secure Com- 
munication 



Finally, wc will discuss very briefly another fascinating concept, namely, syn- 
chronization of chaos, which is attracting considerable interest among chaos 
researchers in recent times. The possibility of two or more chaotic systems oscil- 
lating in a coherent and synchronized way is not an obvious one. One of the main 
features often associated with the definition of chaotic behaviour is the sensitive 
dependence on initial conditions. Then one may conclude that synchronization 
is not feasible in chaotic systems because it is not possible in real systems either 
to reproduce exactly identical initial conditions or exact specification of system 
parameters of two similar systems. We may be able to build "nearly" identical 
systems, but there is always an inevitable technological mismatch and noise, im- 
peding the exact reproduction of all the parameters and initial conditions. Thus, 
even an infinitesimal deviation in any one of the parameters will eventually result 
in the divergence of nearby starting orbits. 

In this connection, the recent suggestion of Pecora and Carrol [20] that it is 
possible to synchronize even chaotic systems by introducing appropriate coupling 
between them has revolutionized our understanding. They have shown that if 
a reproduced part (response or subsystem) of a chaotic system responds to a 
chaotic signal from the full (drive) system, under some conditions the signals 
in the response part would converge to the corresponding signals in the drive 
system. Such a possibility is known as the synchronization of chaotic systems. 
This idea has been further generalized by cascading [20] the reproduced parts or 
subsystems. Here the considered chaotic system may be divided into two sub- 
systems, each of which will synchronize with the full system when driven by the 
appropriate chaotic signal. By arranging one of these synchronized subsystems 
to drive the other subsystem, it is possible to produce synchronized signals in 
the subsystems for every signal in the full drive system. 

Another innovative method which has been developed recently for chaos 
synchronization is the method of one-way coupling between two identical chaotic 
systems [21,22] . In this configuration the response chaotic system variables follow 
identically the drive chaotic system variables for appropriate one-way coupling 
strength. Moreover, the behaviour of the response system is dependent on the 
behaviour of drive system, but the drive system is not influenced by the response 
system's behaviour. 

Thus by reproducing all the signals at the receiver under the influence of a 
single chaotic signal from the drive, the synchronized chaotic systems provide a 
rich mechanism for signal design and generation, with potential applications to 
communications and signal processing. Further, one can use the pair of cascad- 
ing driven subsystems or the identical response system with one-way coupling 
term as attractor recognition devices for chaotic attractors. Hence, a single in- 
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put signal will allow recognition and reconstruction of the driving system and/or 
rejection of non-matching device systems. Another interesting area of applica- 
tion would be cryptography. For example, the non-driven part of the chaotic 
system produces signals, which are typically broad-band, noise-like and difficult 
to predict, and so they can be ideally utilized in various contexts for masking 
information bearing waveforms. They can also further be used as modulating 
waveforms in spread spectrum techniques. 

9.1 Pecora and Ccirroll method 

In this method, Pecora and Carroll have considered an n-dimensional system 
governed by the state equation of the form 

^ = f[X{t)], X = {xuX2,---,Xnf. (53) 

They divide the system into two parts in an arbitrary way as X = {xD^ xr)^ . 
The D part is referred to as the driving subsystem variables and the R part as 
the response subsystem variables. Then Eq. (53) can be rewritten as 

xd = 9{xd,xr), {m— dimensional) (54a) 

iR = h{x£),XBX {k— dimensional) (54&) 

where Xd = {Xi,X2,- ■ ■ ,Xm)'^ , Xr = {Xm+l,Xm+2,- ■ ■ Xn)'^ , 9 = [/l(X),---, 

fm{X)f, h = [fm+i{X), • • • , fn{X)f, and m + fc = n. 

Pecora and Carroll suggested building an identical copy of the response sub- 
system with variables x'^ and drive it with the xjj variables coming from the 
original system. In such a model, we have the following compound system of 
equations: 

Xd = g{xD,xii), {m— dimensional)— drive (55a) 

Xr = h{xD,xii), {k— dimensional)— drive (556) 

x'r = h{xD,x'R), {k— dimensional)— response. (55c) 

Under the right conditions, as time elapses the x'j^ variables will converge asymp- 
totically to the Xr variables and continue to remain in step with the instanta- 
neous values of XR{t). Here the drive or master system controls the response or 
slave system through the xd component. The other component is allowed 
to have different initial conditions from that of xr. The generalized schematic 
set-up for this type of synchronization is shown in Fig. 24. 

One can easily see that if all the Lyapunov exponents of the response sys- 
tem consisting of Eq. (55c) are less than zero, then after the decay of the initial 
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transients, x'j^{t) will be exactly in step with xji{t). The Lyapunov exponents 
for Eq. (55c) are called the conditional Lyapunov exponents (CLE^s). Also 
the response system with negative conditional Lyapunov exponents is called a 
stable subsystem or a stable response system. If at least one of the conditional 
Lyapunov exponents is positive, then synchronization will not occur between the 
x'jf and xji variables. If the x'j^ and x^ subsystems are under perfect synchro- 
nization, then the difference between the x'j^ and xr variables, x|j = — xr, 
will tend to zero asymptotically, that is, limt^oo x*j^{t) 0. Then we have 

x*ji{t) = x'jj{f) - XR{t) = h{xD,x'R) - h{xD,XR). (56) 

If the subsystem is linear, we have 

xUt) = Ax*R{t), (57) 

where A is a (A:x A:) constant matrix. Let the eigenvalues of A be (Ai, A2, • • • , Afe). 
The real parts of these eigenvalues are, by definition, the CL-E's we seek. 

If all of the CL-E's are negative, then lim^^oc x*j^{t) = 0, and the subsystems 
will synchronize; if there is a positive CLE, the subsystems will grow further 
apart as t — > 00 and they will never synchronize. Interestingly, an intermediate 
case occurs if one or more of the CLE^s are zero but none is positive; in this 
case, as t ^ 00, the trajectories of the subsystems will be separated by a fixed 
distance depending upon the initial conditions. 

If the response system is nonlinear in nature, then the conditional Lyapunov 
exponents are not so easily determined and we must resort to numerical proce- 
dures to calculate them. But if the response system is linear (for example, a 
linear circuit with passive elements), it is trivial to calculate the CLE in certain 
cases. The above methodology has been successfully applied to obtain chaos 
synchronization in many important nonlinear systems including the Lorenz sys- 
tem, the Rossler system, the hysteretic circuit, Chua circuit, the driven Chua 
circuit, the DVP oscillator, phase-locked loops, and so on. For applications to 
various systems and applications to secure communications we again refer to [6]. 



10 Conclusions 



In this set of lectiires, we have tried to show how even seemingly simple looking 
nonlinear oscillators can give rise to a very rich and baffling variety of complex 
dynamics involving periodic and chaotic motions. Combined analytical, numer- 
ical and analog simulation studies can help to unravel such complex dynamics 
in a systematic manner. Recent developments also show clearly that chaotic 
dynamical behaviour can be controlled algorithmically by minimal efforts so as 
to avoid any harmful effects of chaos on the structure of the system. Further- 
more, identically evolving chaotic systems can be synchronized in a surprising 
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way through appropriate couphng and the concept can be profitably used for 
spread spectrum communications. Thus the study of even such simple nonlin- 
ear systems can bring out an enormous amount of new knowledge on dynamical 
systems. 

Naturally, it will be of great interest to consider assemblies and arrays of 
such oscillators with suitable couplings. Such studies can lead to interesting 
spatio-temporal patterns and chaos as well as to new structures of practical 
importance. The associated mathematical analysis and physical realizations will 
be areas of intense future investigations. 
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Figure Captions 



Fig. 1(a) (i) Solution curve x{t) of Eq. (2) for u)q = 1.0 (ii) Phase portrait in 

the (x — x) plane corresponding to (i). 

Fig. 1(b) (i) Solution curve x(t) of Eq. (5) for a = 0.5, loq = 1.0. (ii) Phase 
Portrait of (i). 

Fig. 1(c) (i) Solution curve x{t) of Eq.(6) for a = 0.5, / = 1.0, Wq = 1.0, and 
u) = 1.0. (ii) Phase portrait of (i) (transients excluded). 

Fig. 2 Resonance curve of the nonlinear oscillator (8) for (3 > 0. 

Fig. 3 Some typical attractors of a nonlinear oscillator, (a) Trajectories ending 
at a stable focus, (b) Stable limit - cycle attractor of period T and angular 
frequency oj = 2'k/T. (c) A quasipcriodic attractor of the Duffing- van der Pol 
oscillator, (d) Chaotic attractor of the Murali-Lakshmanan-Chua circuit 

Fig.4 A forced LCR circuit (8). Here f{t) = Fssinw^i- 

Fig.5 (a) Wave form v of circuit of Fig. 4. (b) Trajectory plot in the {v — il) 
plane, (c) Trajectory plot in the {v — zl) plane of the experimental circuit of 
Fig. 4. 

Fig. 6 Analog simulation circuit of the driven Duffing oscillator (double-well 
potential). 

Fig. 7 A typical three-segment piecewise-linear characteristic curve of a nonlinear 
resistor (N). 

Fig. 8 Circuit diagram of the simplest dissipative nonautonomous circuit. 

Fig. 9 Various types of steady-state attractors of the Duffing oscillator as a 
function of the forcing strength /. 

Fig. 10 Bifurcation diagrams of the Duffing oscillator in the (/ — x) plane com- 
puted through numerical simulation. 

Fig. 11 (a): (i) Schematic diagram of a reverse bubble, (b): (i) Schematic 
diagram of a bubble, (a): (ii) Bifurcation diagram in the range 0.43 — 0.48 of 
the parameter /. (b): (ii) Bifurcation diagram in the range 0.53 < / < 0.62 

Fig. 12 The (/ — lo) phase diagram indicating regions of different steady states 
exhibited by the double- well potential oscillator fora = 0.1,Wo = — 0.5,/3 = 0.5 
in Eq.(25) with f{t) = fcostot. (vertical lines): LO symmetric; (right slanting 
lines): LO unsymmetric; shaded area: cross- well chaotic motion; cross-thatched 
V areas: SO occurs outside F-shaped regions. Adapted from Ref.[ll]. 
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Fig. 13 Various types of steady state attractors corresponding to Fig. 12. 
Fig. 14 Experimental results of the DufRng oscillator circuit (Fig. 6). (i) Phase 
portrait in the {x — x) plane (ii) Lower trace: solution curve x{t), upper trace: 
forcing. 

Fig. 15 Bifurcation diagram of the BVP Eq.(17) in the ranges (a) Ai G (0, 1.8), 
(b) Ai e (0.6,0.8), and (c) Ai G (1.0,1.3). 

Fig. 16 The complete devil's staircase for Eq.(17) at Ai = 0.88 with a = 0.8, 
b = 0.7, c = 0.09, and u = 1.0 

Fig. 17 MLC circuit (28): (a) Numerical simulation results: Phase portraits 
in (x - y) plane: (i) / = 0.065, Period-1 limit cycle; (ii) / = 0.08, Period-2 

limit cycle; (iii) / = 0.091, Pcriod-4 limit cycle; (iv) / = 0.1, One-band chaos; 
(v) / = 0.15, Doubleband chaotic attractor. (b) Experimental results: Phase 
portrait v versus Vs{= Rsih)' (i) F = 0.0365Kms, Period-1 limit cycle; (ii) 
F = 0.0549V;ms, Pcriod-2 Umit cycle; (iii) F = O.OGWrms, Period-4 limit cycle; 
(iv) F = 0.072Wrms, One-band chaos; (v) F = 0.0107Vi.ms; Doubleband chaotic 
attractor. 

Fig. 18 Codimension one bifurcations (see Table 2). 
Fig. 19 Codimension two bifurcations (see Table 3). 

Fig. 20 Regions of LO attractor: computer simulation and theoretical stability 
limits: a = 0.1 (Ref.[ll]). 

Fig. 21 (a) Square lattice of poles foimd in typical solution of free undamped 
DufRng oscillator (a = 0, / = 0). (b) Singularity distribution in the complex 
f-domain of a solution to the damped oscillator. 

Fig. 22 Typical singularity distribution in complex t-domain for a solution of 

the forced Duffing oscillator. Accentuated poles (heavy dots) indicate commonly 
observed 'tunnel' structures (adapted from Ref.[13]). 

Fig. 23 Evolution of the variable x of the BVP oscillator, (a) Limit cycle motion 
without control; (b) Attraction to the fixed point {xs,ys) = (—1.001,-0.376) 
under the control for e = 0.1. Conversion of (c) chaotic motion to (d) a period 
-T limit cycle and (e) period- 2T limit cycle. Aif is the value of Ai associated 
with the desired solution. 

Fig. 24 Schematic representation of cascading synchronization. Drive and re- 
sponse #1 (left slanting lines): Variables used to drive the response #1. (right 
slanting lines): The replica part of drive system used in response, (blank): Drive 
variables which are fed directly into the response ^1. Response =f/=2. Variables 
which are fed directly into the response #2 from response #1. (left slanting 
lines): The replica part of the drive system. 
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Table 1: Summary of bifurcation phenomena of Eq. (25) with a = 0.5, 
(jg = -1.0, /3 = 1.0 and lj = 1.0. 



Value of / 


Nature of solution 


Attractor in the phase 
space in the {x,x) plane 
(simulation results) 
Numerical analog 
Fig. 9 Fig. 14 


/ = o 


Damped oscillation to the 
stable focus at x < 


(a) 


< / < 0.34 


Period-T oscillation 


(b) 


0.34 < / < 0.358 


Period- 2T oscillation 


(c) 


0.358 < / < 0.362 


Period-4T oscillation 


(d) 


0.362 < / < 0.37 


Further period-doubling 




0.37 < / < 0.4 


One-band chaos 


(e) 


0.4 </ = 0.42 


Double-band chaos 


(f) 


0.42 < / 


Chaos, windows, reverse 
period-doublings, chaos, 
boundary crisis, etc. (See 
bifurcation diagram of 
Fig. 10) 


(g), (h) 
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Table 2: Codimension one (one-parameter family of) bifurcations 



Type of bifurcation 


Model equation 


Nature of bifurcation 


Flows 
1. 


Saddle node 


(i) i = /i — 

(ii) X = fi + x'^ 
a: e 7^^ At e 7^l 


Fig. 18a 
Fig. 18b 


2. 


Transcritical 


X = fj, — x^ 


Fig. 18c 


3. 


Pitchfork 

(i) Supercritical 

(ii) Subcritical 


X = fix — x^ 
X = fix + x'^ 


Fig. 18d 
Fig. 18e 


3. 


Hopf 

(i) Supercritical 


r = r{fi — r^) 
e = 1 
=x'^ + 


Fig. 18f 




(ii) Subcritical 


f = r{fi + r^) 
6 = 1 


Fig. 18g 


Maps 
4. 


Flip 

(period doubling) 


X G [0,1), < // < 4 


Fig. 181i 



Table 3: Codimension two (two-parameter family of) bifurcations 



Type of bifurcation 


Model equation 


Nature of bifurcation 


1. 


Cusp catastrophe 


a; = Ai -|- \2X — x^ 
x e TV-, (Ai,A2) e TV 


Fig. 19a 


2. 


Degenerate Hopf 


f — r(Ai -|- A2r^ -|- ar^) 
e = u) + 0{r'^) 


Fig. 19b 


3. 


Takens -Bogdanov 


x = y 


Fig. 19c 






y = \ix + X2y - x^ - x^y 
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